Octupolar ordering of classical kagome antiferromagnets in two and three dimensions 
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Classical Heisenberg antiferromagnets on two-dimensional kagome and three-dimensional hyperk- 
agome lattices are investigated by Monte Carlo simulations. For both models the symmetry-breaking 
states at low temperatures are described by non-zero octupole moments or third-rank spin tensor 
order parameters. In the case of the two-dimensional kagome antiferromagnet, a sharp crossover 
into a coplanar state takes place at Tk w 0.004J, which we attribute to proliferation of fractional 
vortices. The three-dimensional model exhibits a first-order transition at T c ~ 0.002J into a phase 
with critical spin correlations, which possesses a long-range order of octupole moments. 

PACS numbers: 75.10.Hk, 75.10.-b, 75.50.Ee, 75.40.Mg 



I. INTRODUCTION 

A two-dimensional (2D) network of corner-sharing tri- 
angles known as the kagome lattice, Fig. [T] is a proto- 
type of geometrical frustration. The nearest-neighbor 
Heisenberg antiferromagnet on such a lattice has an in- 
finite number of spin configurations minimizing the ex- 
change energy. Both quantum^ and classical 2 -i£ii£i£i 7 .i& 
spin models on the kagome lattice have attracted sig- 
nificant theoretical interest in the past. Realizations 
of the kagome lattice topology among magnetic solids 
were initially rather scarce with the prime example being 
SrCrg_ a; Ga4+ a ;Oi9i^ii^ In the last few years a significant 
number of new magnetic compounds that are believed 
to be related to the kagome lattice antiferromagnet have 
been synthesized and studied^ 2 ^ 14 ! 15 ' 16 ! 17 . 18 ! 19 . 20 ! 21 
Often these materials suffer from substitutional disorder, 
are affected by small structural deviations from the ideal 
kagome network, or have extra interactions, which lift the 
magnetic degeneracy. Nevertheless, recent neutron scat- 
tering experiments on powder samples of large- S kagome 
materials, Yo.sCao.sBaC^Oy— (S = 3/2) and deutero- 




nium jarosite (S = 5/2), have demonstrated remark- 
able similarity between the measured diffuse intensities 
and the Monte Carlo results for the classical models 
Motivated by these two seemingly good realizations of 
the classical kagome antiferromagnet, we reinvestigate in 
the present work the finite-temperature properties of this 
model. In particular, we consider the angular depen- 
dence of magnetic correlations, which can be measured 
in neutron-diffraction experiments on single crystals. 

A second source of motivation is provided by the re- 
cent discovery of a 3D array of corner-sharing triangles 
m a spin- 1/2 Mott insulator N& 4 lr 3 8 H Due to simi- 
larity with its 2D counterpart, this lattice structure has 
been coined a hyperkagome lattice. A network of tri- 
angles with similar topology is also know to exist in 
gadolinium gallium garnet GdaGasOia- 2 ^^ (5 = 7/2), 
whose enigmatic behavior attracted a lot of theoretical 
efforts ; 26 ' 27 ! 28 Though the magnetic properties of both 
systems may be quite distant from those of the nearest- 
neighbor classical model, in the former case due to quan- 
tum effects, in the latter material due to strong dipo- 



FIG. 1: (Color online) Section of the kagome lattice with spins 
in the fully ordered v3 x \/3 structure. 



lar interactions, it is still important to understand the 
properties of the classical antiferromagnet as a start- 
ing reference point. Furthermore, recent Monte Carlo 
simulations^ 9 - have found evidence for an interesting low- 
temperature phase transition for the hyperkagome anti- 
ferromagnet. 

The classical ground states of the Heisenberg kagome 
lattice antiferromagnet are derived from the block repre- 
sentation of the spin Hamiltonian 



n 



j 



^(S! + S 2 + S 3 



A 



const . (1) 



The energy is minimized by any spin configuration, which 
has Sa = for every triangular plaquette. This classi- 
cal constraint is satisfied for infinitely many configura- 
tions including planar and nonplanar states. Chalker et 
al£ have argued that coplanar spin states are selected 
by thermal fluctuations via the order by disorder effect. 
They have also related an asymptotic selection of the spin 
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plane with development of the nematic orde r 30 i 31 in spin 
chiralities defined as 

k = ^=(Si x S 2 + S 2 x S 3 + S 3 x Si) (2) 

for each triangular plaquette. The published Monte 
Carlo data seem to confirm this prediction 2 ^ and the cor- 
responding point of view prevails now in the literature on 
frustrated magnets. 

Below we present arguments that such a description is 
incomplete and the low-temperature state of the classical 
kagome antiferromagnet should be described by a third- 
rank tensor or octupolar order parameter. Such a pro- 
posal was first put forward a long time ago^^ though 
no numerical results were presented to substantiate this 
idea. The difference between the broken symmetries for 
the two types of order parameters is important for topo- 
logical classification of point defects in the kagome anti- 
ferromagnet. Topologically stable defects or vortices play 
a significant role in low temperature transformations of 
2D geometrically frustrated magnets and may lead to 
topological phase transitions 32 and/or to an unconven- 
tional spin-glass behavior^ The analytic consideration 
is supported in the following by extensive Monte Carlo 
simulations. 

The paper is organized as follows. In Sec. II the pos- 
sible tensor order parameters are considered for magnet- 
ically disordered spin systems and presence of fractional 
vortices is emphasized in the case of the classical kagome 
antiferromagnet. Section III is devoted to Monte Carlo 
results for the 2D kagome antiferromagnet. In particular, 
the specific heat exhibits a sharp kink, which signifies for- 
mation of the coplanar spin state. In Sec. IV we investi- 
gate the behavior of the 3D kagome antiferromagnet and 
find that it shows a fluctuation-driven first-order transi- 
tion. The low-temperature phase possesses no long-range 
antifcrromagnetic correlations and is described instead 
by an octupolar order parameter. 

II. SPIN TENSOR ORDER PARAMETERS 
A. Symmetry analysis 

Let us first recall the arguments for the order by dis- 
order effect in the kagome antiferromagnet Similar 
discussion for the hyperkagome antiferromagnet is post- 
poned until Sec. IV. The classical constraint Sa = 
for one triangular plaquette is satisfied by a 120° spin 
structure. Once the orientation of the three sublattices is 
chosen for the first plaquette, all lowest energy coplanar 
configurations can be identified with the ground states 
of the three-state Potts antiferromagnet or, equivalcntly, 
with coloring all sites of the kagome lattice into three 
colors, such that no two neighboring sites have the same 
color. There are 1.13471^ such states with N being the 
number of lattice sitesi^ 3 . The manifold of coplanar states 
contains two simple periodic structures: the q = state, 



in which all triangles pointing up (A) or down (y) are 
in the same state, and the \/3 x y/3 structure, which 
is shown in Fig. Q] and is described by the wave-vector 
Q = (47r/3o, 0), a being the lattice constant. The q = 
state has the same chiralities for up and down triangles, 
while in the y/3 x V3 structure chiralities alternate be- 
tween A and y plaquettes. 

Nonplanar states are constructed from planar config- 
urations by identifying various closed or open two-color 
lines. In the y/3 x V3 state, these are represented by 
hexagonal loops, see Fig. [TJ Nearest neighbors off such a 
line are necessarily spins of the third color. Spins on the 
line can be continuously rotated about the direction de- 
termined by the third sublattice. The obtained spin fold 
(also called weather-vane mode) retains the 120° spin 
orientation and costs, therefore, no energy. Rotation by 
7r returns spins back into a single plane creating a new 
coplanar state. 

The harmonic analysis^ indicates that coplanar states 
are selected at low temperatures because they have the 
largest number of soft excitations. The harmonic excita- 
tion spectra arc identical for all coplanar configurations. 
Hence, selection of a specific translational pattern, if any, 
occurs due to weaker nonlinear effects. Therefore, there 
should be a range of temperatures where selection of the 
spin plane is not accompanied by a wave- vector selection. 
The heat capacity in this regime is equal to C = 11/12 
per spin3 The spin plane is specified by its normal ([2]), 
which for a general chirality-disordered coplanar state se- 
lects a line without direction. The corresponding order 
parameter is a second-rank traceless tensor; 30 ! 31 

p 

where summation extends over all triangular plaquettes. 
A simpler form of the nematic order parameter can be 
constructed as a sum of on-site quadrupole moments: 

i 

The two order parameters [Eqs. © and describe the 
same type of broken symmetry and it is only a matter of 
convenience to choose one of them. 

This is not, however, the end of the story. The copla- 
nar states break, in addition, the spin-rotational sym- 
metry inside the plane: At large distances spins do 
not follow any specific translational pattern but still are 
chosen from the initial sublattice triad. The ground 
states of the XY kagome antiferromagnet with planar 
spins Sj = (cos 8 'j, sin 8 'j) have a long-range order in 
wj = exp(3i#j)4 Generalization to Heisenberg spins is 
given by an on-site octupole moment expressed as a sym- 
metric third-rank tensor: 

ifT = S?S? 57 - \ SfSfH - \ Sf6 ai - \ S7S afi (5) 

with vanishing trace over any pair of indexes. The uni- 
form long-range order of such octupoles is described by 
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non-zero values of 

rpc 



(6) 



where (• • •) denotes thermodynamic averaging. The ten- 
sor T a ^ has in total seven independent components, as 
follows from its symmetry and tracelessness. Note, that a 
similar duality in the choice between the two tensor order 
parameters exists for liquid crystals consisting of bent- 
core molecules.— ^ A complete characterization of the 
orientational order in such systems requires definition of a 
third-rank tensor order parameter in addition to the more 
familiar nematic tensor. Different forms of third-rank 
spin tensors have been discussed in the literature i 5 ' 8 ' 36 i 37 
For classical spins all of them are equivalent to Eq. |J5]), 
the latter form being more convenient for numerical sim- 
ulations. 

The order parameters T 01 ® 1 and Q a @ transform accord- 
ing to different irreducible representations of the rota- 
tion group SO (3) corresponding to the angular momenta 
I = 3 and 1 = 2, respectively. This does not mean, how- 
ever, that octupole and quadrupolc moments cannot co- 
exist below the critical point. The two order parameters 
are coupled by a rotationally invariant term in the free- 
energy functional 



(7) 



Due to time-reversal symmetry, the octupole moment has 
zero average value in the nematic phase. In contrast, an 
instability driven by T"^ 1 also induces a nonvanishing 
quadrupolar tensor unless (^"-^T 13 ^) = 8 a p. Since the 
ensemble of coplanar state is described by nonzero val- 
ues of both tensors, the primary order parameter for the 
kagome antiferromagnet is the octupole moment ^ . The 
quadrupolc moment in this case is only a secondary order 
parameter. In the nomenclature of phase transition the- 
ory, the low-T state of the kagome antiferromagnet can 
be called 'improper spin nematic.' Numerical data in sup- 
port of the above conclusion are presented in Sec. IIIB. 

Possible symmetries of a uniform octupolar state are 
deduced by minimizing the Landau free-energy func- 
tional: 

_|_ v jiaffj rpaS/J, rpfiSu ijr^iu ^ 

For v > 0, the stable phase has D3 symmetry, whereas for 
v < 0, it is invariant under tetrahedral point group i 34 i 35 
The low-temperature state of the kagome antiferromag- 
net is, naturally, identified with the D3-symmetric "tri- 
atic" state: D3 spin rotations permute three spin sublat- 
tices and transform one translationally disordered copla- 
nar state into another one from the same ensemble. The 
spin tensor for the triatic state can be parametrized as 

T a ^ cx (l a fV - rm^m 7 - m^m 7 - m a m l3 P) (9) 

with two orthogonal unit vectors 1 and m lying in the spin 
plane. The "tetrahedratic" phase (v < 0) is not realized 



in the present spin model. Note, that the interaction 
term ([7]) vanishes in the tetrahedratic state. Hence sta- 
bilization of the triatic phase can be ascribed to a strong 
coupling between the octupolar and the nematic order 
parameters in the kagome antiferromagnet. 

For 2D Heisenberg antiferromagnets the true long- 
range order is impossible at any finite temperature. The 
above discussion applies in this situation to a symme- 
try of spin correlations at short distances. For dis- 
tances larger than the correlation length r > £ the or- 
der parameters (J3])— d6j) vanish and one has to consider 
instead the generalized susceptibilities xv — N(V 2 )/T 
with V = Q a @ and T a ^ 7 . The lattice-averaged squares 
of the two spin tensors, which are directly measured in 
Monte Carlo simulations, are 



and 



(T 



^^[((s.-s,) 3 )-^^^,; 

i ,3 



(10) 



(11) 



At zero temperature, in the fully ordered triatic phase 
(the ground state of the three-state Potts model), the 
above expressions yield the following limiting values 
((Q«/3) 2 ) = 1/6 and ((T Q ^ 7 ) 2 ) = 1/4 . 



B. Topological analysis 

A state with ordered octupole moments also exhibits 
nontrivial topological properties^ The order-parameter 
space of the kagome antiferromagnet is obtained as a 
coset space R = SO(3)/Z?3. It has the non-Abelian ho- 
motopy group tt\ (R) = D 6 . The symmetry of the octupo- 
lar phase allows, therefore, the presence of stable topolog- 
ical defects. It is instructive to compare the topological 
properties of the kagome antiferromagnet with those of 
its weakly frustrated counterpart, the triangular lattice 
antiferromagnet. The latter spin system is conventionally 
ordered into the 120° spin state, which breaks completely 
the rotational symmetry. The long-range translational 
order leaves no discrete symmetries of the spin struc- 
ture. As a result, the degeneracy space of the order pa- 
rameter is R = SO (3) and point defects are ^-vortices: 
7Ti[SO(3)] = A possible realization of ^-vortex is 

shown in Fig. [2^ and corresponds to a pattern with all 
spins lying in a common plane and performing a 2-7T ro- 
tation around the vortex center. Equivalently, it can be 
represented as rotation of equilateral triangles with dis- 
tinguishable vertices. 

For a coplanar state of the kagome antiferromag- 
net, spin folds permute three spin sublattices and the 
slow spatial variations of the magnetic order parame- 
ter are represented by a texture of equilateral triangles 
with equivalent vertices. An elementary point defect 
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(a) (b) 

FIG. 2: (Color online) Stable topological defects in non- 
collinear Heisenberg antiferromagnets. The three-sublattice 
120° spin structure is locally represented by orientations of 
equilateral triangles. Two examples of the spin texture cor- 
respond to (a) Z2-vortex in the triangular lattice antiferro- 
magnet with the winding number n — +1 and (b) fractional 
vortex in the kagome antiferromagnet with n = +1/3. 

is, consequently, a fractional vortex with ±2-7r/3 rota- 
tion around the vortex core, see Fig. These 1/3- 
vortices are known to reduce substantially the Kosterlitz- 
Thouless transition temperature in the XY kagome 
antiferromagnet i 39 ' 40 

In Heisenberg magnets topological defects may have a 
more complicated nonplanar structure. Simple hydrody- 
namic arguments suggest that the nonplanar Z2-vortex in 
the triangular antiferromagnet has a lower energy com- 
pared to the planar vortex in Fig. ^.M* The same line 
of arguments applies also to defects in the kagome an- 
tiferromagnet. In addition, the hydrodynamic energy of 
1/3- vortices is further reduced by a factor of 1/9 due to 
the smaller phase winding. As a result, entropic genera- 
tion of fractional vortices in the kagome antiferromagnet 
starts at significantly lower temperatures than a similar 
effect for ^-vortices in the triangular antiferromagnet. 
A possible role of the non-Abelian topological defects in 
the kagome antiferromagnet was brought to attention in 
Ref. d and will be further discussed in Sec. IIIB. 



III. KAGOME ANTIFERROMAGNET 

A. Monte Carlo algorithm 

The published Monte Carlo data for the nearest- 
neighbor Heisenberg antiferromagnet on the kagome lat- 
tice were performed on relatively small clusters of N — 
3L 2 spins with L < 24 (N < 1728)M>! Numerical re- 
sults in the present work have been obtained for a sub- 
stantially wider range of lattices with L = 12- 72. The 
standard Metropolis algorithm has been adopted. A site 
on a periodic cluster is randomly picked up and a new 
orientation of spin is chosen. The new direction is ac- 
cepted according to the Metropolis rejection scheme. To 
increase acceptance rate a maximum variation AS Z = T 
on a z-component of spin in the local coordinate frame is 



imposed at low temperatures. In this way the acceptance 
rate stays close to 50% in the whole temperature range. 
A sweep over the lattice in which on average every spin 
is attempted to move corresponds to one Monte Carlo 
(MC) step. 

For the Heisenberg kagome antiferromagnet the slow- 
ing down develops into a serious problem at low temper- 
atures T/J < 0.01. The autocorrelation time can be fur- 
ther reduced by using the microcanonical over-relaxation 
algorithm^ Generally, for Heisenberg models the over- 
relaxation move consists in successive rotations of spins 
around their respective local field by an arbitrary angle 
such that the total energy remains unchanged. The sim- 
plest and most efficient realization corresponds to a it 
rotation, i.e., flipping a spin to the most distant direc- 
tion from the initial one4£ Implementation of such a spin 
move requires neither generation of random numbers nor 
calculation of trigonometric functions, which saves signif- 
icantly operation time. Lattice scans can be performed 
with random or sequential selection of spins. We find that 
collective motion of spins is generated more efficiently 
in the latter procedure yielding faster decorrelation. Fi- 
nally, one hybrid MC step consists of one canonical MC 
step followed by a few microcanonical nonrandom up- 
dates. Such deterministic reshuffling of spins is essential 
for reducing autocorrelation times at low temperatures. 
Typically we use between 3-10 over-relaxation updates 
per one MC step depending on cluster size. 

Each finite cluster was initiated with a random spin 
configuration and gradually cooled to the lowest tem- 
perature T/J = 10 -4 . At every temperature 5 x 10 4 
hybrid MC steps were allowed for equilibration which 
were followed by measurements (~ 5 x 10 5 ) performed in 
intervals of five hybrid MC steps. In addition, all mea- 
sured quantities have been averaged over 20-50 cooling 
runs, starting from different random configurations. This 
further helps to overcome a freezing problem and also 
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FIG. 3: (Color online) Temperature dependence of the specific 
heat for a kagome lattice cluster with L = 36. The horizontal 
arrow denotes the value C/N = y|, the two vertical arrows 
indicate boundaries between three different regimes. 
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FIG. 4: (Color online) Finite-size behavior of the specific heat 
in the low temperature region. 



provides an unbiased estimate for the statistical errors. 
Unless otherwise specified, the error bars do not exceed 
the symbol sizes. Special checks have been performed to 
verify that the hybrid MC algorithm works efficiently in 
the relevant temperature range when instead of gradual 
cooling we start from either a random spin configuration 
or the ordered \/3 x V3 structure. Full thermal equili- 
bration of the q = ground state was achieved only for 
T/J> 0.002, which is still significantly better than in 
the previous studies £*L 



B. Macroscopic properties 

Let us begin with the heat capacity, which has been 
computed from fluctuations of the internal energy C — 
{{E 2 ) — (E) 2 )/T 2 . The temperature dependence of 
C(T) on a linear-logarithmic scale is shown in Fig. [3] 
for a cluster with L = 36. One can clearly distin- 
guish three different regimes for the specific heat with 
the two crossover points indicated by vertical arrows. 
The high-temperature regime T/J > 0.1 corresponds 
to a paramagnetic phase with only weak correlations 
between neighboring spins. In the intermediate regime 
0.005 <T/J< 0.1 the internal energy reaches its classi- 
cal minimum value £/JV = - J up to a small contribu- 
tion from thermal fluctuations. Spins on triangular pla- 
quettes become strongly correlated and satisfy approxi- 
mately the constraint condition Sa = 0. This regime is 
commonly known as a classical spin-liquid or a coopera- 
tive paramagnet42, The specific heat in the cooperative 
paramagnetic state remains close to C/N — 1, which re- 
flects absence of soft modes in the excitation spectrum. 

Selection of smooth, locally coplanar spin configura- 
tions takes place at T/J < 0.005 as indicated by a re- 
duced specific heat. The probability distribution peaks 
in the vicinity of coplanar ground states, which have one 
zero (anharmonic) mode for every hexagon. The limiting 
value C/N\t^o coincides quite accurately with 11/12 = 



0.916 . . . predicted by the mode counting analysis^ 

The enlarged low-temperature part of C(T) is shown 
in Fig. 3] for several cluster sizes. Two features are note- 
worthy. First, the crossover between a planar spin state 
and a cooperative paramagnet corresponds to a rather 
sharp kink in C{T) at around T k /J = 0.004 ± 0.0005. 
At T < Tk the specific heat grows linearly with tempera- 
ture, which can be accounted for by interaction between 
the spin waves. Second, the specific heat exhibits a pecu- 
liar finite-size behavior in the vicinity of the kink point, 
showing a rounded peak on small clusters, which disap- 
pears for L > 30 with no significant finite-size corrections 
afterward. 

The temperature dependence of the mean-square of 
the octupole moment is shown in Fig. [5^,. Large clus- 
ters exhibit a clear enhancement of the order parameter 
below Tk/J ~ 0.004, which coincides with a kink posi- 
tion in the specific heat. At low temperatures ((T Q ^ 7 ) 2 ) 
approaches 1/4, which is the limiting value for the fully 
ordered coplanar phase. The octupolar susceptibility 
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FIG. 5: (Color online) Temperature dependences of (a) the 
octupolar order parameter and (b) the corresponding suscep- 
tibility for different cluster sizes. 
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FIG. 6: (Color online) Temperature dependence of the three 
order parameters normalized to their respective values in the 
fully ordered state. Numerical data are for a cluster with 
L = 36. 



is presented in Fig. [SJd. For each cluster there is an in- 
flection point T£ below which the correlation length £r 
becomes of the order of the linear lattice size £r ~ L 
and the susceptibility begins to exhibit finite-size effects. 
The lattice independent part of xt diverges as T — ► 
signaling a long-range ordered state at T = 0. The fast 
increase in xt{T) at low temperatures is consistent with 
a typical divergence xt{T) — AT n exp(B/T) found from 
the nonlinear sigma model mapping^ though no specific 
predictions for A, B, and n exists for the kagome antifer- 
romagnet. Note, that a rapid crossover in the behavior 
of xt(T) takes place in the vicinity of T/J ~ 0.005. 

The behaviors of the two order parameters T Q ^ 7 and 
Q a @ is compared in Fig.[6l The octupole moment shows a 
faster growth with decreasing temperature, which would 
correspond to a larger exponent if a second-order tran- 
sition is assigned to Tfc. Fig. [5] illustrates our previous 
conclusion that the octupolar order parameter drives the 
low-temperature transformation in the kagome antiferro- 
magnet. 

Finally, we investigate the elastic properties of the 
kagome antiferromagnet by computing the temperature 
dependence of the spin stiffness. The spin stiffness p s 
is defined as the second derivative of the free-energy 
with respect to weak nonuniform twist of spins performed 
about a certain direction a in spin space: 



AF = - 
2 



d 2 rp s [VQ a {v)f 



(13) 



Substituting 9f — S8(e-Ti) for the twist angle and taking 
the limit 80 — > we obtain the following expression after 
proper symmetrization and normalization per unit area: 



V3 ri 

~2N\'3 



(E) 



£!/[^( Si xS j r(e-5, 



(ij) 



)}■ 

(14) 



where (E) is the internal energy and e is an arbitrary 
unit vector in the lattice plane. 

Numerical results for p s (T) are presented in Fig. [7] for 
three system sizes. The spin stiffness vanishes at temper- 
atures above T/J ~ 0.005. This further supports identi- 
fication of the intermediate phase at 0.005 <T/J< 0.1 
with the cooperative paramagnet, which has well devel- 
oped local spin correlations but exhibits zero response to 
long-wave-length perturbations. The spin rigidity starts 
to increase at T/J < 0.004 simultaneously with the devel- 
opment of short-range octupolar correlations. Finite-size 
scaling in the low temperature regime still yields p s — 
as it should be for a 2D Heisenberg spin system. 

The origin of the sharp crossover in various proper- 
ties of the kagome antiferromagnet at T ~ Tk ~ 0.004 J 
deserves special attention. The possible phase transi- 
tion in 2D continuous non-Abelian models driven by 
nontrivial topological defects has been discussed in the 
context of two different physical applications. The first 
group of works motivated by investigation of liquid crys- 
tals has studied the RP 2 model in 2D, which in the 
spin language corresponds to a model of three-component 
spins on a square lattice coupled with ferro-biquadratic 
exchangei 45 i 46 i 47 ' 48 The order parameter space is the 
projective sphere RP 2 — S 2 /Z2 with the first homo- 
topy group n\{RP 2 ) = The topological defects in 
this context are called disclinations. Independently, the 
role of topological defects was emphasized for 2D non- 
collinear Heisenberg antiferromagnets . 32 ' 49 '^^ 1 -'^'^ The 
order parameter space is SO(3) = S 3 /Z2 in this case and 
the fundamental group is the same 7Ti[SO(3)] = Z2. For 
both types of models a straightforward generalization of 
the Kosterlitz-Thouless scenario suggests that topologi- 
cally stable Z2-vortices are bound in pairs for T < T v 
and become free in the high-temperature phase d 2 ' 45 

Kawamura and Miyashita 32 investigated numerically 
the Heisenberg antiferromagnet on a triangular lattice 
and found evidence for the vortex unbinding transition at 
T v ~ 0.3 J. The heat capacity exhibits a weak maximum 
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FIG. 7: (Color online) Spin stiffness of the kagome antiferro- 
magnet for different cluster sizes. 
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FIG. 8: (Color online) Temperature dependence of the anti- 
ferromagnetic order parameter for different cluster sizes. Hor- 
izontal arrows indicate corresponding values in the ground 
state of the three-state Potts model. 



in the vicinity of T v . The main difference with the stan- 
dard Kosterlitz-Thouless transition in planar spin sys- 
tems is that the correlation length remains finite both 
above and below T v . This leads to a small finite density of 
free defects in the low temperature phase. The low- and 
the high-temperature phases are still distinguished by an 
asymptotic behavior of the vorticity on a large closed 
contour: the vorticity function changes from the perime- 
ter law at T < T v to the area law at T > T V M£&& The 
renormalization group analysis becomes, however, noto- 
riously difficult since in this case it must include simul- 
taneously spin-waves and Z2-vortices. The precise form 
of a singularity in the thermodynamic potential at such 
a topological transition remains unknown up to now^i 

The topological properties of the kagome antiferro- 
magnet suggest a natural interpretation of the observed 
crossover in terms of unbinding of fractional vortices. It 
may also provide another example of topological tran- 
sition in 2D Hcisenberg antiferromagnets. The kink 
anomaly in C{T) is consistent with a cusp-type sin- 
gularity in the specific heat found at the topological 
transition for the RP 2 model^ and for the triangular 
antiferromagnet *22, The behavior of the spin stiffness also 
agrees with the defect unbinding scenario. Similar to the 
Kosterlitz-Thouless transition, free topological defects 
are responsible for vanishing p s — above the crossover 
point, whereas a much slower decrease in p s with the sys- 
tem size at low temperatures is determined by spin-wave 
excitations. Further numerical studies, which directly 
measure the density of fractional vortices and the corre- 
sponding vorticity function, are necessary to clarify the 
above conjecture of topological transition in the kagome 
antiferromagnet . 



C. Spin correlations 

The high-temperature series expansion for the kagome 
antiferromagnet 3 finds that the maximum in the 
momentum-dependent susceptibility corresponds to the 
spin structure. A similar conclusion has been 
made by Huse and Rutenberg^ from a different perspec- 
tive: spin correlations of the three-state Potts model 
arc dominated by the staggered component at the wave- 
vector of the V3 x v3 structure with a power-law decay 
~ r -4/3 a t i on g distances. Such a purely entropic effect 
derived from the mapping to the two-component height 
model is related to the fact that the 'flat' y/3 struc- 
ture maximizes the number of flippable loops. An en- 
hancement of the antiferromagnetic correlations at low 
temperatures was also seen in the Monte Carlo simula- 
tions of the Heisenberg modeli^ 

We have investigated spin correlations in the kagome 
antiferromagnet by using the following form of the stag- 
gered magnetization: 



m 2 AF 



-y 

N 2 ^ 



jQ(Ri — Rj) 



(15) 



where index I numbers spins in the unit cell, i, j and Rjj 
denote position of the unit cell on the triangular Bra- 
vais lattice, and Q = (47r/3, 0) is the wave-vector of the 
\/3 x \/3 structure. The lattice constant, which is twice 
the intersite spacing, is chosen as the length unit. The 
normalization factor gives m 2 ^ = 1 in the fully ordered 
structure. The temperature dependence of the antiferro- 
magnetic order parameter for several lattice sizes is pre- 
sented in Fig. [H Using the loop-flip algorithm^ we have 
also measured the antiferromagnetic order parameter for 
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FIG. 9: (Color online) Distance dependence of spin correlators 
normalized to the correlations of the v3 x v3 structure. Cor- 
relators are measured along chain directions on the kagome i 
lattice of linear size L = 48. Distances are given in units of 
the lattice constant. 
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the three-state Potts antiferromagnet at T — 0. The cor- 
responding values are shown in Fig. [5] by the horizontal 
arrows. 

Spin correlations may also yield sharp signatures in the 
static magnetic structure factor S(q); 

5 (q) = ^E^- s ^ eiq(r ^ rj) . ( 16 ) 

i,3 

where is a spin position in the 2D plane. The energy in- 
tegrated neutron scattering cross-section is proportional 
to the instantaneous spin correlator (|16p and provides an 
experimental tool for its measurement. Previous numeri- 
cal works on the kagome antiferromagnet have considered 
only the powder averaged structure factorial 

In the whole range of temperatures the antiferromag- 
netic order parameter has much smaller values than the 
octupole or quadrupole moments, see also Fig. [6l The 
enhancement of m\ F observed at T/J < 0.01 on small 
systems^ is significantly suppressed for larger clusters. 
For every lattice size L there is a characteristic tempera- 
ture below which the amplitude of the -\/3 x y/3 correla- 
tions becomes larger than the corresponding correlations 
in the ground state of the three-state Potts antiferro- 
magnet. This happens because the 7r-folds dressed with 
short wave-length fluctuations acquire a finite linear ten- 
sion, whereas loops of all lengths are flipped with equal 
probability in the three-state Potts model. Still, the en- 
hancement relative to the three-state Potts model is pro- 
gressively shifted to lower temperatures with increasing 
cluster size and in the end might be a finite-size effect. 

Explicit comparison of spin correlations in the Heisen- 
berg antiferromagnet and the three-state Potts model is 
presented in Fig. [9] for a lattice with L — 48. A simi- 
lar plot is given in Fig. 5 of Ref. 0. Once the difference 
in the length scales is taken into account, the previously 
studied cluster corresponds to L = 12 in our notations. 
At T/J = 0.005, the spin correlations on the L = 12 lat- 
tice approach or even exceed the amplitudes of the Potts 
model* 4 - which is consistent with the enhancement of the 
staggered magnetization shown in Fig. [SJ On the other 
hand, the correlations for the L = 48 cluster at the same 
T fall consistently below the Potts model amplitudes, see 
Fig. [9] There is, therefore, a fundamental difference be- 
tween the low-temperature behaviors of the spin tensor 
order parameters and ttiaf- While the octupole suscep- 
tibility exhibits an exponential growth as T — > 0, which 
signifies a finite value of the order parameter at T = 0, 
the staggered susceptibility xaf has a much weaker in- 
crease, which for available system sizes does not exceed 
corresponding values of the three-state Potts antiferro- 
magnet. 

Over the last several years it has been established that 
certain 2D and 3D classical spin models governed by local 
constraints exhibit critical behavior with power-law de- 
cay of spin correlations at long distances i 4 ' 54 i 55 ' 56 i 57 ' 58 i 59 
A widely discussed consequence in application to the py- 
rochlore antiferromagnet is a "bow-tie" shape of the mag- 
netic structure factor with pinch points in certain high- 



symmetry directions i 58 ' 59 i 60 Expansion around the large- 
N limit for a 0(N)-symmetric spin model on the kagome 
lattice also yields a 1/r 2 decay law for spin correlations 
with similar bow-tie features in S(q)~ It is, therefore, 
interesting to compare magnetic intensities in the three 
different temperature regimes of the Heisenberg kagome 
antiferromagnet (Fig. [3]) with the above predictions. 

The magnetic structure factor has been calculated for 
the kagome lattice cluster with L = 24 at four different 
temperatures. The obtained results are shown in Fig. |TD] 
as intensity plots on a square slice in the reciprocal space 
with the extension from to 87r in both directions. The 
highest temperature T/J — 0.5 corresponds to the para- 
magnetic state, where thermal fluctuations are strong 
and the magnetic intensity has only a broad structure in 
the momentum space. In the cooperative paramagnetic 
state at T/J = 0.1 the pinch points develop between tri- 
angular shaped regions of strong intensity. In the vicinity 
of the pinch point S(q) has a nonanalytic form due to 
long-distance spin-spin correlations with dipolar-like an- 
gular anisotropyi 55 i 58 ' 59 Upon cooling to a lower temper- 
ature T/J = 0.02, the intensity is redistributed in favor 
of centers of the triangular regions, which correspond to 
Q' = 2Q = (87r/3, 0) and equivalent wave- vectors. These 
are not true Bragg peaks: as the spin correlations fall off 
as r~ 2 in this regime^ the peak intensity grows loga- 
rithmically with the system size. Note, that peaks at 
the antifcrromagnetic wave-vector Q and equivalent po- 
sitions are absent. 

Finally, when the coplanar correlations start to develop 
at T/J = 0.005, the narrow necks loose significantly in 
intensity, while new satellite peaks at q — Q becomes 
noticeable. The intensity maps for T/J — 0.1, and 0.02 
most closely resemble the analytic result for S(q) from 
the large-TV expansion^ However, the magnetic inten- 
sity exhibits more structure once the order by disorder 
effect selects coplanar states. The development of ex- 
tra diffuse peaks in the low-temperature regime is also 
observed in the powder averaged structure facto n 7 ' 19 ? 21 
Note, that the algebraic decay of spin correlations dis- 
cussed above occurs at distances smaller than the cor- 
relation length. The inverse correlation length provides 
a natural width for all the nonanalytic features in the 
magnetic structure factor 5(q). 



IV. HYPERKAGOME ANTIFERROMAGNET 

The best known example of 3D geometrically frus- 
trated lattice is a network of corner-sharing tetrahedra of 
the pyrochlore lattice, see, e.g., Ref. [6ll The experimen- 
tal realizations include numerous magnetic pyrochlore 
and spinel compounds. The only example of a 3D lat- 
tice of corner-sharing triangles was so far provided by 
gadolinium gallium garnet GdaGasOi^ . 23124 The recent 
experiment has found another interesting example of a 
3D triangular network of magnetic ions in Na^^Os- 22 
The corresponding lattice structure can be obtained by 
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FIG. 10: (Color online) Magnetic structure factor of the kagome antiferromagnet in different temperature regimes. Region in 
the momentum space corresponds to < q x , y < 8tt. Temperatures from left to right and from top to bottom are T/J — 0.5, 
0.1, 0.02, and 0.005. The intensity scales differ by a factor of 2 between the rows. 



1/4 depletion of the pyrochlore lattice such that only 
three out of the four vertices of each tetrahedron are oc- 
cupied by magnetic ions, see Fig. 1111 By analogy with 
the kagome lattice, this structure is called the hyperk- 
agome lattice. It contains 12 spins in the standard cubic 
unit cell. The garnet lattice has, in contrast, 24 spins 
in its cubic unit cell and consists of two interpenetrating 
hyperkagome sublattices. The positions of the atoms in 
each garnet sublattice are different from the 1 / 4-depleted 
pyrochlore structure but topology of the two networks re- 
mains the same. The local symmetry on magnetic sites 
is somewhat higher in the garnet structure containing a 
two-fold rotational axis joining centers of adjacent trian- 
gles. 

In this Section we consider a classical Heiscnberg an- 
tiferromagnet on the hyperkagome lattice. The limited 
applicability of such a model to the above two magnetic 
materials has been mentioned before. Still this model is 
quite interesting on its own, in particular, to contrast it 
with the more familiar 2D system. The arguments given 
before for the infinite degeneracy of the classical ground 
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FIG. 11: (Color online) The hyperkagome lattice as a 1/4- 
depleted pyrochlore structure. 



states of the kagome antiferromagnet (Sec. I and II) are 
fully applicable to the 3D model and will not be repeated 
here. The coplanar configurations play again an impor- 
tant role at low temperatures. The harmonic analysis 
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FIG. 12: (Color online) Low-temperature specific heat of the 
classical hyperkagome antiferromagnet. 




FIG. 13: (Color online) Mean-square of the octupolar order 
parameter versus temperature for the hyperkagome antiferro- 
magnet. 



finds N/3 zero-energy modes for the hyperkagome anti- 
ferromagnet out from the total 27V modes for N classical 
spinsi 27 ' 29 This yields the same limiting value for the spe- 
cific heat C/N\ T ^ = H/12 as for the 2D model. Soft 
modes reside on closed loops, which pass through at least 
ten triangles. 27 ' 62 

A submanifold of coplanar configurations is again 
mapped to the ground states of the three-state Potts an- 
tiferromagnet on the hyperkagome lattice. The precise 
number of such states or, equivalently, the ways of color- 
ing sites with three colors is not known for the 3D lattice. 
A local gauge representation of the classical constraint 
with the Maxwellian-type action predicts the dipolar 
form of the spin-spin correlation function! 29 ' 57 Perform- 
ing loop-algorithm simulations in the ground state en- 
semble of the Potts antiferromagnet we indeed found a 
very fast decay of spin correlations consistent with an 
r~ 3 dependence. 

The general mode counting arguments suggest 
the presence of the order by disorder effect for 
three-component spins on a lattice of corner-sharing 
triangles. 61 Though the early numerical work found no 
evidence of such an effect^ the recent Monte Carlo study 
detected a first-order transition at T/J < 0.004.— The 
authors argued in favor of a nematic state below the 
transition. As we have shown in Sec. II the full sym- 
metry breaking pattern in the submanifold of the copla- 
nar states is described by the octupolar order parameter. 
The first-order nature of the transition is not surpris- 
ing in this respect as the renormalization group analysis 
of the Landau free-energy functional {H finds no stable 
fixed-point solution below four dimensions^ 

To verify the formation of the octupolar ordering in 
the hyperkagome antiferromagnet we have performed the 
classical Monte Carlo simulations on periodic clusters 
with N — 12L 3 spins and L = 3-6. The equilibra- 
tion problem at low temperatures in three dimensions 
becomes much more severe compared to the 2D model. 



The published data 2 ^ show, for example, a strong hys- 
teresis for L > 4, which clearly indicates that large clus- 
ters fall out of equilibrium. We resort, therefore, to the 
exchange Monte Carlo algorithm 63 in a conjunction with 
the hybrid updates described before. We use between 
30 to 50 replicas depending on cluster size distributed in 
the temperature interval 0.001 < T/J < 0.02. All repli- 
cas are initiated with random spin configurations, which 
are equilibrated for 10 6 exchange MC steps. These are 
followed by measurements for 5 x 10 5 MC steps. Finally, 
the results have been averaged over 20 independent runs 
to determine statistical errors. 

The numerical data for the specific heat obtained on 
two smallest clusters are presented in Fig.[T2J At temper- 
atures below the anomaly, the specific heat approaches 
the value C/N — 11/12 predicted for the coplanar state 
by the mode counting analysis. The fast growth of the 
peak height is consistent with C/N cx L 3 scaling ex- 
pected for the first-order transition^ 

The temperature dependence of the mean square of 
the octupolar order parameter is shown in Fig. 1131 The 
transition temperature for each cluster can be estimated 
at the mid-height of the jump. The order parameter ex- 
hibits extremely weak finite size dependence at low tem- 
peratures, which signifies development of the true long- 
range octupolar ordering in this 3D spin model. Still, a 
certain amount of disorder in the low-temperature phase 
is evidenced by deviations of ((T Q/37 ) 2 ) from the limiting 
value of 1/4 for the fully ordered octupolar structure. 
A simple scaling of the transition temperature with the 
cluster size yields T c /J = 0.002 ± 0.0003 for the tran- 
sition temperature. In order to locate more accurately 
the point of first-order transition, one needs to simulate 
bigger lattices, which appears to be a difficult task even 
for the exchange Monte Carlo algorithm. 

In order to check a tendency to a long-range magnetic 
ordering below T c at a priori unknown wave vector we 



11 



1.5 



0, 





1 i 1 


/ "A 




/ A 
I \ 


- 


J \ * 


V f 


// Xl^ 


- 


1 








i — T/J 


= 0.002 


// T/J 


= 0.2 


/ / 




— 


i 



FIG. 14: (Color online) Powder averaged magnetic structure 
factor for the hyperkagome antiferromagnet. Wave-vectors 
are given in reciprocal lattice units. 

have calculated the angular-averaged structure factor: 

s( „4 ?(s ,. s ,igM. (17) 

The neutron diffraction experiments on polycrystallinc 
samples allow to directly measure S(q). The results 
for the L = 5 cluster are presented in Fig. [TJ] The 
static structure factor S(q) below T c exhibits broad dif- 
fuse peaks with no sharp features, which could point at 
a long-range order. We have found only a tiny change in 
S(q) across the first-order transition. Variation of tem- 
perature by 2 orders of magnitude from T/J = 0.2 to 
0.002 leads only to a small increase ~ 15% in the main 
peak intensity. A similar form of the diffuse structure 
factor at high temperatures was previously obtained in 
the numerical work on gadolinium gallium garnet^ The 
nearest-neighbor classical Heisenberg antiferromagnet on 
the hyperkagome lattice provides a unique example of the 
long-range ordered state of magnetic octupoles with al- 
gebraically decaying spin-spin correlations. 

Returning back to Na^^Os, the recent experiment 22 
has demonstrated the absence of magnetic ordering down 
to 2 K, which is significantly smaller than the scale of 
antiferromagnetic interactions deduced from the Curie- 
Weiss constant dew ~ 650 K. Such a behavior is con- 
sistent with strong geometrical frustration found for the 
nearest- neighbor hyperkagome antiferromagnet 1 26 i 29 i 62 It 
would be interesting to compare the neutron diffraction 
data in the spin-liquid state of Na^^Os with the above 
results for S(q). This should allow one, in particular, to 
verify relevance of the nearest-neighbor Heisenberg model 
to the real material. 



V. CONCLUSIONS 

We have clarified in the present work that selection 
of the coplanar states in the kagome antiferromagnet 



is properly described by the development of the oc- 
tupolar (third-rank spin tensor) order parameter. Our 
Monte Carlo simulations yield Tk ~ 0.004J for the on- 
set of coplanar ordering, which is lower than the previ- 
ous estimates Furthermore, we suggest that Tk may 
correspond to a topological transition, which consists in 
unbinding of fractional vortices. Presence of these topo- 
logically stable point defects follows from the nontrivial 
degeneracy space of the octupolar order parameter. The 
antiferromagnetic \/3xy/3 correlations are also enhanced 
at low temperatures, though the corresponding correla- 
tion length remains significantly shorter than the charac- 
teristic length scale for octupolar correlations. Our MC 
data for big lattices with L > 36 do not confirm the pre- 
viously made suggestion- that the asymptotic \/3 x \/3 
ordering develops in the limit T — > 0. Precise numeri- 
cal study of the lowest temperature region T/J < 10~ 3 
definitely requires simulations on clusters with L > 72, 
which is impossible without new Monte Carlo algorithms 
(but see also the last paragraph). 

The classical hyperkagome antiferromagnet provides a 
unique example of the spin model with the long-range 
octupolar ordering. In the broken-symmetry state be- 
low a fluctuation induced first-order transition spin cor- 
relations remain critical with a power-law decay r~ 3 at 
long distances! 29 ' 57 It is interesting to study the phase 
diagram in magnetic field for this model, in particular, 
in relation to the experimental diagram of Gd3Ga5 0i2 
which exhibits a field-induced ordered phase i 23 ' 24 ' 25 

Let us finish with a few comments on the possible re- 
lation of the topological transition to the spin-glass like 
behavior observed in some jarositesi 12 ' 13 Similar to the 
previous numerical studies we find a very strong tendency 
to spin freezing once the coplanar configurations are sta- 
bilized at T < T v . It is the hybrid Monte Carlo algo- 
rithm adopted in the present work, which allows one to 
equilibrate large spin systems at low temperatures. The 
spin reshuffling dynamics of the microcanonical sweeps 
has no simple analog in real magnetic materials. The 
kagome antiferromagnet will be, therefore, stuck in one 
of the many degenerate coplanar states with frozen struc- 
ture of chirally ordered domains. Rotation of spins along 
spin folds is suppressed due to development of free-energy 
barriers of the entropic origin^i 5 -^ Above the topological 
transition, when thermally excited vortices destroy the 
common plane and remove the entropic barriers, the sin- 
gle spin-flip dynamics becomes effective again. The above 
ideas for the unconventional spin-glass transition were pi- 
oneered by Ritchey and co-workersi 5 ^ The results of our 
work put emphasis on a hidden topological transition 
behind a spin-glass freezing in geometrically frustrated 
magnets. Though the corresponding temperature scale 
comprises only a small fraction of 9qw for the Heisen- 
berg antiferromagnet, intrinsic XY anisotropies in real 
materials can significantly enhance T v by transforming it 
to the Kosterlitz-Thouless transition for 1/3- vortices. 
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